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Abstract. Some of the contributions of Chandrasekhar to the field of magnetohydrodynamics are 
highlighted. Particular emphasis is placed on the Chandrasekhar-Kendall functions that allow a 
decomposition of a vector field into right- and left-handed contributions. Magnetic energy spectra 
of both contributions are shown for a new set of helically forced simulations at resolutions higher 
than what has been available so far. For a forcing function with positive helicity, these simulations 
show a forward cascade of the right-handed contributions to the magnetic field and nonlocal inverse 
transfer for the left-handed contributions. The speed of inverse transfer is shown to decrease with 
increasing value of the magnetic Reynolds number. 
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1. Introduction 

During the 1950s, Chandrasekhar was deeply immersed in studying problems in hydro- 
dynamics and especially magnetohydrodynamics (MHD), which is the study of flows of 
electrically conducting media such as liquid metals and ionized gases or plasmas. Be- 
tween 1950 and 1960, roughly 2/3 of his papers were related to hydrodynamic and mag- 
netohydrodynamics stability, turbulence, and aspects of dynamo theory. Chandra's scien- 
tific activities during that period are reviewed in detail by Parker 1 1], who draws a close 
connection between Chandra's interests in Heisenberg's theory of turbulence [2] and the 
then discovered polarization of starlight in the galaxy ii. This led to a collaboration 
between Chandrasekhar and Fermi [5], where they estimated the galactic field strength to 
be 6-7 /xG, which is close to the commonly accepted value today L6]. Parker also suggests 
that observations of magnetic fields in the Sun |7] and other stars 081] may have contributed 
to Chandra's interest in magnetohydrodynamic stability and perhaps dynamo theory. 
Chandrasekhar's name is not normally connected with dynamo theory. Indeed, at the 

time there was still great disbelief in dynamo theory, given that Cowling |@] and others, 

I — I 

including Chandra himself [10], proved dynamos impossible under the assumption of ax- 
isymmetry of the magnetic field. After all, Chandra's work in that field was concentrated 
to the period just before Herzenberg lHHl produced the first existence proof of dynamos. 
Therefore, Chandrasekhar kept his mind open for alternative explanations. In 1956 this 
led him to the question [12] "Can the time of decay of 14,000 years, in the absence of 
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internal motions, be prolonged to 500,000 years (say) by velocities of reasonable magni- 
tude and patterns?", which he answered tentatively with "Yes". He discussed this finding 
in connection with possible solutions to the problem of the geomagnetic field and even 
the field observed in sunspots. Now we know that the physically correct and relevant 
understanding, at least for the Sun, came from Parker's work \ ]3\ of 1955. However the 
special significance of Parker's early work was perhaps not yet fully appreciated when 
Chandrasekhar lll2ll in 1956 merely noted that "the possibility, nevertheless, of dynamo 
action has been explored, intensively, in recent years by BuUard & Gellman 1 14], Parker 
m, and by Elsasser Elil]." 

Chandrasekhar's most cited publication is his book on Hydrodynamic and hydromag- 
netic stability which, in 2010, received about 140 citations on ADS. The citation 
rate has been increasing and was about 70 citations some 10 years earlier He studied 
the stability of virtually any imaginable combination of systems in slab and spherical ge- 
ometries, with and without rotation, with and without magnetic fields, as well as systems 
with self-gravity. In astrophysics, hydrodynamic and hydromagnetic instabilities are im- 
portant sources of turbulence. Unlike atmospheric and wind tunnel turbulence, which 
is typically driven through boundary layer instabilities, turbulence in astrophysical set- 
tings tends to be driven by instabilities within the volume. Important examples include 
Rayleigh-Benard convection as well as the magneto-rotational instability. Chandrasekhar 
has contributed to both types of instabilities. It is now well accepted that at large enough 
Reynolds numbers, the resulting flows tend to become turbulent and are invariably subject 
to yet another instability: the dynamo instability. Of particular significance is the case of 
large-scale dynamo action. All known mechanisms produce helical large-scale magnetic 
fields that tend to be nearly force-free. This leads to yet another of Chandra's favorite 
topics at the time, namely the study of such fields in spherical geometry in terms of what 
is now known as Chandrasekhar-Kendall functions 1 18]. In the present paper we discuss 
the significance of these functions for understanding the dynamics of magnetic fields ex- 
hibiting inverse transfer behavior from small to large length scales. This mechanism is 
significant not only in stars , accretion discs ]21], and galaxies ll22ll . but possibly 



even on the scale of the Universe in connection with primordial magnetic fields 0231 12411 . 



2. Chandrasekhar-Kendall functions and application to Cartesian domains 



When it was realized that astrophysical magnetic fields might often be force-free 112511 . i.e. 
(V X B) X B = 0, Chandrasekhar began studying the mathematical properties of such 
fields both in cylindrical and spherical coordinates 126J . Such fields are eigenfunctions of 
the curl operator, i.e., 

V X hn= knhn, (1) 

where the fc„ are eigenvalues and /i„ eigenfunctions. These functions play important 
roles as basis functions for solenoidal vector fields as well as for decomposing them into 



right- and left-handed components. This led to an important joint paper with Kendall II18I1 . 
whose supervisor (V C A Ferraro) remarked in a footnote of their paper that "The results 
in this paper were derived independently by the two authors; and they agreed to write it as 
one". Their work gained in importance since the 1970s, when force-free magnetic fields 
were observed in the solar corona [27] and produced in laboratory plasmas [28]. 

Chandrasekhar-Kendall or CK functions continue to be used in solar physics, for exam- 
ple in describing the dynamics of coronal loops [29]. These functions have subsequently 
been exploited also in other geometries, including cylindrical geometries [30], where they 
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are also referred to as Lundquist fields iJsT] . In the plasma physics community, these func- 
tions are sometimes also referred to as Chandrasekhar-Kendall-Furth or CKF functions 



1 32I1 . In the following, however, we will focus on the Cartesian domains with periodic 
boundaries i33ll . where such functions are also known as Beltrami waves. 

In periodic Cartesian domains it is convenient to adopt CK functions in Fourier space. 
For example the magnetic vector potential can be written in Fourier space (indicated by 
subscript fc) as 

Ak{t) = j A{x,t)e'''-^d?x, (2) 

and it can be expanded in terms of eigenfunctions of the curl operator with positive and 
negative eigenvalues, as well as a longitudinal part that vanishes upon taking the curl to 
compute the magnetic field B = 'W x A. Thus, we write i24l33ll3411 



with 



Au{t) - a+{t)h+ + a^{t)h-^ + al{t)hl (3) 



ikxh^^±kh^, k^\k\, (4) 



where are the eigenfunctions with positive and negative eigenvalues, respectively. 
These functions are orthonormal, i.e., 

(K* ■h^) = {K* ■hl) = {h^* ■hl)=0 (5) 

and 

(K* ■K) = {hu* ■K) = {hV -^,1) = !, (6) 

where asterisks denote the complex conjugate, and angular brackets denote volume av- 
erages. The longitudinal part a^h^l is parallel to k and vanishes after taking the curl (or 
crossing with ifc) to calculate the magnetic field, i.e., we have 

Bk = k {a+h+ - a^h^) . (7) 

Magnetic helicity and magnetic energy spectra are then given by 

Ak-Bl = k[ia+f -ia^r-], (8) 



B^-Bl = kmatr + {a~r]. (9) 

This implies evidently 

Bk-Bl>k\Ak-Bl\, (10) 

which is also known as the realizability condition fsS*], that is often quoted in terms of 
shell-integrated spectra of magnetic energy and magnetic helicity, Mk and Hk, respec- 
tively, i.e., 

Mk>k\Hk\/2fio. (11) 

Here, /io is the vacuum permeability and the 2/io factor comes from the normalization 
of the magnetic energy spectrum, J Mkdk = (B^/2/io), while the magnetic helicity 
spectrum is normalized such that J Hk dk ^ {A ■ B) . 

In order to monitor the chiral dynamics of a system, it is convenient to plot energy 
spectra of the contributions from right- and left-handedly polarized parts. Fortunately, 
this can be done without ever performing the explicit decomposition Q in terms of a^. 
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Figure 1. Visualizations of Bx/Bcq on the periphery of the domain at three times 
when the small-scale magnetic field has already saturated. Here, Boq ~ y/fiopo Wrms 
is the equipartition field strength where kinetic and magnetic energy densities are com- 
parable, and po is the mean density. Note that the maximum field strength is about 
twice Ben ■ 



3. Application to the inverse cascade 

Eigenfunctions of the curl operator are well suited as function basis for hydrodynamic tur- 
bulence |,33i1 . More recently they have also been used in magnetohydrodynamics [24], and 



in dynamo theory 03411 . where flows allow a conversion of kinetic into magnetic energy. 

Dynamos are possible under completely isotropic conditions without helicity. In that 
case we often speak about small-scale dynamos; see Ref. |36] for a review. However, in 
the presence of helicity, there is an inverse cascade of magnetic helicity that leads to the 
gradual build-up of magnetic fields at the largest scale of the system \ l^,^2B]. This has 
been demonstrated in real space with turbulence driven by an ABC-flow forcing function 
(3T\\ and, more clearly, with turbulence driven by a forcing function proportional to CK 
functions that are (5-correlated in time (3^. In the latter case, the inverse cascade pro- 
cess is clearly demonstrated in terms of CK-decomposed spectra, M^, as was done in 



Ref. 113411 . They used a relatively big scale separation ratio, i.e., their forcing wavenum- 
ber kf was large compared with the smallest wavenumber of the domain, ki (they used 
kf/ki = 27). This implies, however, that the magnetic Reynolds number, 

-Rm = Urms/??fcf, (12) 

is smaller than if the scale separation was less. It is therefore useful to increase the numer- 
ical resolution as much as possible. In this paper we present new simulations of helically 
driven isotropic turbulence in a periodic domain at a resolution of 512'^ meshpoints and 
scale separation ratios of kf/ki — 15 and 30. The resulting magnetic Reynolds number 
reaches values up to « 57 and the magnetic Prandtl number is unity. The simulations 
have been performed with the PENCIL COD^E which is a modular high-order code (sixth 
order in space and third-order in time, by default) for solving a large range of partial 
differential equations. 

In Figure [T] we show the x component of the magnetic field at the periphery of the 
domain at two early times in the simulation. Note, that the magnetic field appears to 
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Figure 2. Spectra of magnetic energy of the negatively polarized contributions, M^T , 
at earlier (left) and later (right) times. The scale separation ratio is fcf/fci = 15. The 
range of time tis given in units of the turnover time, t — 1/ "Umis fcf . At small wavenum- 
bers, the A/jT spectrum is proportional to fc*, while to the right of fcf /fci — 15 there is 
a short range with a. k^^ spectrum. 



organize itself gradually into larger scale patches. This is more clearly seen in energy 
spectra; see Figure |2] where we show the spectral magnetic energy at different times. At 
early times the magnetic energy resides almost entirely at the dissipative scale. How- 
ever, as time goes on, a magnetic field component at intermediate scale with wavenumber 
km ~ fcf/2 becomes discernible. This corresponds to the fastest growing eigenmode of 
a purely helical large-scale dynamo ['sT]. Note also that near the largest scale of the sys- 
tem, the magnetic field has a power spectrum comparable with k*, which corresponds to 
white noise in the magnetic vector potential and obeys causality i39ll . To the right of the 
injection wavenumber, kf, there is a short range of a k^^ spectrum in . 

At later times, magnetic energy is transported further to smaller wavenumbers; see the 
right-hand panel of Figure |2] where the peak of the spectrum at k = k^ travels to smaller 
wavenumbers. The effect of this is also evident in real space in that larger scale patches 
of magnetic field with the same orientation occur In Figure |3] we show the x component 
of the magnetic field at the time when a large-scale magnetic field has emerged, but it has 
still not reached the scale of the computational domain. Magnetic energy spectra of that 
time are shown in the upper panel of Figure |4] Note, however, that the magnetic energy 
does show a clearly pronounced secondary peak that is now at wavenumber fcm/ ki « 2. 

The process that produces magnetic fields at large length scales shown in the present 
paper can loosely be regarded as an inverse cascade of magnetic helicity and thereby 
magnetic field. However, it can also be treated in terms of mean-field theory, in which 
case it corresponds to the so-called a effect, whose physics was already identified by 
Parker lll3ll . This dynamo is quite different from the laminar dynamos that were studied 
at the time by Bullard and Gellman [[Till , although there remained substantial uncertainty 
about the convergence of the numerical solutions obtained back then. 

One of the surprising results of the last 10 years is the fact that in a closed or periodic 
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Figure 3. Visualizations of Bx/Bcq on the periphery of the domain at a later time 
when a large-scale magnetic field is beginning to emerge. This time corresponds to the 
intermediate time shown in the second panel of Figure|2]where the peak of the spectrum 
has reached the wavenumber fcm/fci ~ 3, corresponding to about 3 large-scale patches 
across the domain. The scale separation ratio is fcf/fci — 15, so there are about 15 
eddies across the domain. 



domain, the large-scale magnetic field saturates only on a resistive time scale fss"]. This 
is understood as a consequence of the magnetic helicity equation, which limits the time 
scale for magnetic helicity evolution to the resistive time 140(1 . and lowers the nonlinear a 
effect in an i?ni-dependent fashion . According to this picture, if we associate, as 

was done in Ref. 13811. the peak wavenumber fcm of the inversely cascading field with the 
a effect via = a/2rit, where i-jt is the turbulent magnetic diffusivity, we should expect 
the speed of inverse transfer to slow down as well. When this was first proposed [34], the 
results were not fully conclusive and of course in conflict with the usual understanding 



that the inverse cascade should operate on a dynamical time scale II45I1 . Our present results 



allow us now a fresh look at this problem. As a continues measure of fcm we define 

^m" = / Mk dk / [ Mk dk , (13) 



where we choose n = 3, which was found to give the best approximation to the actual 
location of the secondary peak in Mk- This is probably related to the sharp rise of Mk 
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Figure 4. Upper panel: spectra of magnetic energy, M^, and rescaled magnetic helic- 
ity, ±kHk/2iio. Lower panel: spectra of magnetic energy of positively and negatively 
polarized parts, . 



proportional to k'^ for fci < fc < fcm. Following earlier work ll34ll . the evolution of kj^it) 
is well described by the fit formula 

k~^ = Q;trav(i - tsat), (14) 

where the parameter atrav characterizes the speed at which fcm travels toward smaller 
wavenumbers. This parameter is not directly related to the a effect, although atrav also 
has the dimension of velocity. The result is shown in Figure |5] where we plot fcj"^ versus 
time for different values of i?i„ and fcf /fci, and we also show the dependence of the slope 
atrav on _Rm. In agreement with earlier work [jsH, atrav is clearly dependent and 

— 1/2 

decreases with increasing i?,„ roughly proportional to Rm ■ The reason for this seems 
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Figure 5. Left panel: time dependence of the peak wavenumber for scale separation 
ratios of 15 (dashed) and 30 (solid lines) at Rm of 12, 27, and around 57 (increasing 
line thickness). Right panel: _Rni dependence of the cascade speed for scale separation 

— 1/2 

ratios of 15 (open symbols) and 30 filled symbols. The straight lines give the 
(dotted) and (dashed) dependences. 



clear: the bump in the spectrum is a manifestation of an a effect causing nonlocal spectral 
transfer to smaller k, but it is not an inverse cascade in the usual sense. An inverse cascade 
may still occur at intermediate wavenumbers fcm ^ fc ^ fcf , where evidence for constant 
(fc-independent) spectral magnetic helicity transfer to larger scale has been found Ii46ll . 

4. Closing the loop: turbulence from the magnetic field itself 

Large-scale dynamos show mean magnetic fields that can either be statistically steady, 
or they can oscillate and propagate in space. In the language of mean-field dynamos, 
the former ones are typically the so-called dynamos while the latter ones are the ail 
dynamos, where the ft effect corresponds to differential rotation or shear amplifying the 
toroidal magnetic field from a poloidal one. A prominent example of such a dynamo is 



that believed to exist in accretion discs, where shear plays a strong role 112111 . This brings 
us to our final example where Chandrasekhar's contributions are still being mentioned in 
modern astrophysics, namely in connection with the magneto-rotational instability. 

Already in 1953, Chandrasekhar studied the hydromagnetic stability of Couette flows, 
e.g. the flow of a liquid metal between two rotating cylinders with angular velocities fti 
and ft2 at radii i?i and R2 > Ri- In view of astrophysical applications to the stability of 
accretion discs, particularly interesting is the case where angular velocity decreases with 
radius, i.e. fti > 5^2, but angular momentum increases, i.e. R^ili < i?2^^2- Under the 
assumption of a small magnetic Prandtl number, where the viscosity is much smaller than 



the magnetic diffusivity, and a small gap width, Chandrasekhar Il47[l found that the flow 
is stable both with and without an axial magnetic field. However, relaxing the assump- 
tion of a small gap width, Velikhov 148.1 found that for an ideal fluid (zero viscosity and 
zero magnetic diffusivity) the flow is unstable when ili > £72. This was the also con- 
firmed by Chandrasekhar |49], but the astrophysical significance was hardly appreciated 
until Balbus & Hawley |50] rediscovered this instability, which is now usually called the 
magneto-rotational instability or MRl. 

The point is that both Velikhov and Chandrasekhar were very much ahead of their 
time, because accretion discs entered the astrophysical literature only in the late 1960s 
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and early 1970s. Nevertheless, it is important to realize that in the meantime the MRI was 
not completely forgotten. In his book, Safronov f5]\ discussed several mechanisms for 
causing turbulence in pwtostellar discs. Among the various mechanisms, he did discuss 
the MRI, but discarded it on the grounds that protostellar discs are cold and poorly ionized, 
so magnetic fields would not be tied to the gas. His work was before the famous paper of 
Shakura & Sunyaev lf52ll . where the physics of accretion discs was applied to the much 
hotter discs around stellar-mass and supermassive black holes. Curiously enough, until 
then the focus has rather been in stability, as is evidenced by a quote from S H Hong 1531] 
who writes "A special form of this problem was first investigated by Chandrasekhar ll47ll . 
and later by Velikhov in i48ll . On the basis of the small gap approximation, they found 
that the effect of a sufficiently strong magnetic field is to inhibit the onset of instability." 

The MRI is now generally believed to be the main agent that drives turbulence in accre- 
tion discs, and this turbulence is able to reinforce the magnetic field by dynamo action via 
a doubly-positive feedback , Sin. These discoveries would not have been possible 



without the help of supercomputers. However, coming back to Chandrasekhar, it is fair to 
say that he always had an eye toward numerical solutions. For example his approach to 
the dynamo problem in 1956, where he found a significant slow-down of the decay fl^l, 
was also entirely based on numerical solutions. It is therefore clear that his combined 
approach using numerical and analytical tools was not only a modem one, but much of 
his work attracts interest still today, as is evidenced by the increasing citation rate of his 
book on Hydrodynamic and hydromagnetic stability Iil7i1 . 
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Allocations Committee. This work was supported in part by the European Research Coun- 
cil under the AstroDyn Research Project No. 227952 and the Swedish Research Council 
Grant No. 621-2007-4064, and the National Science Foundation under Grant No. NSF 
PHY05-51164. 
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